MnO Spin- Wave Dispersion Curves from Neutron Powder Diffraction Data 
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We describe a model-independent approach for the extraction of spin-wave dispersion curves from 
neutron total scattering data. The method utilises a statistical analysis of real-space spin configura- 
tions to calculate spin-dynamical quantities. The RMCProfile implementation of the reverse Monte 
Carlo refinement process is used to generate a large ensemble of supercell spin configurations from 
powder diffraction data. Our analysis of these configurations gives spin-wave dispersion curves that 
agree well with those determined independently using neutron triple-axis spectroscopic techniques. 
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In this Letter, we explore the unexpected possibility 
that neutron powder diffraction might be used to probe 
spin dynamics in magnetic materials. In effect, what 
we are asking is whether there is sufficient information 
preserved in experimentally observed neutron scattering 
functions to allow reconstruction of the spin- wave disper- 
sion. 

It is well known that the intensities of magnetic Bragg 
reflections are determined entirely by one-particle cor- 
relations (the same is true, naturally, for the nuclear 
Bragg reflections). In terms of spin-dynamical informa- 
tion, these describe the average fluctuations in moment 
at each magnetic site. However, any determination of 
spin-wave dispersion relations would require access to the 
two-particle spin correlation functions. For this reason, 
we have focussed our efforts on an analysis of the total 
magnetic neutron scattering patterns, considering both 
elastic (which includes Bragg) and inelastic (included in 
the diffuse) components alike 0, Q . This is because the 
latter cannot be described effectively without recourse 
to the two-particle correlation functions; indeed, all the 
spin-dynamical information is folded into the total scat- 
tering process. What is not clear is the extent to which 
the orientational averaging implicit in powder diffraction 
experiments complicates — or indeed precludes — recovery 
of this information. Our paper aims to address the two 
issues so raised: (i) to what extent are the two-particle 
spin correlations preserved in total neutron scattering 
data, and (ii) how might one extract these correlations in 
practice, and in doing so determine spin-wave dispersion 
relations? The potential impact of a robust methodol- 
ogy based on neutron powder diffraction is significant, 
particularly since spin- wave information could be rapidly 
obtained from e.g. newly-discovered polycrystalline ma- 
terials and/or those materials for which a route to large 
single crystal samples does not exist. 

A similar problem — one that has received more 
widespread attention — is that of extracting phonon dis- 
persion curves from the neutron powder diffraction pat- 
terns of non- magnetic materials 0, 0, Q On the one 



hand, the extraction of spin- waves might be expected to 
prove more tractable, in that the number of spin-wave 
modes in a material is always less than the number of 
phonon modes, and the absolute energies of spin- waves 
tend to be lower than most phonon energies (and hence 
their signature within diffraction data should be clearer) . 
On the other hand, the analysis of diffraction data from 
magnetic materials will be complicated both by the su- 
perposition of nuclear and magnetic scattering contribu- 
tions and by the magnetic form-factor dependence. What 
is clear — and this has emerged in the development of 
methods for extracting phonons from diffraction — is the 
need for model- independent approaches, so that the an- 
swers one obtains are driven by data rather than the 
choice of any particular dynamical model 0, El • 

This paper describes a "proof of principle" investiga- 
tion into the possibility of extracting spin-wave disper- 
sion curves from powder diffraction data. We present 
results for manganese(II) oxide (MnO) that indicate one 
can indeed recover the spin- wave excitation spectrum for 
this system, albeit with some identifiable yet surmount- 
able limitations. Our paper begins by describing the re- 
verse Monte Carlo (RMC) method of refining total neu- 
tron scattering data. We proceed to illustrate how sta- 
tistical analysis of large ensembles of these RMC con- 
figurations can yield quantitative measurements of the 
spin-wave dispersion. Our diffraction-based MnO spin- 
wave dispersion curves are presented and compared with 
those obtained independently using triple-axis inelastic 
neutron spectroscopy. We close with a discussion of some 
necessary developments for the approach to offer a viable 
method of measuring spin- wave dispersion curves in ma- 
terials for which established techniques are impractical. 

For non-magnetic materials, the Fourier transform of 
the observed total scattering data gives the real-space 
atomic pair distribution function, which contains detailed 
information regarding the local structural environments 
present in the sample. The RMC method provides a 
particularly effective means of analysing total scattering 
data, in that it can be used to refine atomistic configura- 
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tions that represent large supercells of the known struc- 
tural unit cell 0, Q . These supercells reflect simultane- 
ously the average periodic structure (as given by the ob- 
served Bragg intensities) together with the collective dis- 
placements from which the diffuse scattering arises. We 
have shown elsewhere that these atomic displacements 
are well described as an instantaneous superposition of 
all phonon modes 0, El • 

Magnetic materials require refinement of both atom- 
istic and spin configurations to account for the observed 
scattering function. We have recently extended the 
RMCProfile implementation of the RMC method to allow 
simultaneous modelling of nuclear and magnetic struc- 
tures ||- The procedure, which is described in detail 
elsewhere [4, 8] , involves minimisation of the "mismatch" 
function 

2 _ 2 , 2 / -1 \ 
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where 
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Here, S(Q) is the observed scattering function as mea- 
sured in each of k data sets (such as from a range of de- 
tector banks, for example), and I(t) is the time-of-flight 
Bragg powder profile as described in more detail else- 
where |7(. The function Xbmc * s minimised by succes- 
sive random moves of two sorts: displacement of atoms 
in the atomistic configuration and re-orientation of indi- 
vidual spin vectors in the associated spin configuration. 
Such moves are continually generated and accepted or 
rejected according to the Monte Carlo algorithm until 
the fit-to-data does not improve further. At this point, 
one has a pair of configurations (one atomistic, one spin) 
that reflects both nuclear and magnetic contributions to 
the Bragg and diffuse scattering patterns. As such, we 
can expect the RMC spin orientations to represent an 
instantaneous superposition of the spin-wave modes: a 
"snapshot" of the spin dynamics. The essence of our 
analysis is to collect large numbers of these spin configu- 
rations, so as to sample many different "snapshots" of the 
spin-dynamical motion. By probing the various spin cor- 
relations that recur, and the frequency with which they 
do so, we will show how one can reconstruct spin- wave 
dispersion curves from such ensembles. 

We have elected to use MnO as a case-study in this 
investigation as its magnetic structure and spin dynam- 
ics have been the focus of a large body of systematic 
investigation || 0, 0, 0, 0] . ft is widely considered 
a benchmark antiferromagnetic material, and is com- 
monly used as a representative of the family of first-row 
transitional-metal oxides in a variety of theoretical stud- 
ies. A paramagnet at room temperature, MnO exhibits 



long-range antiferromagnetic ordering at temperatures 
below Tn = 118 K. Its paramagnetic to antiferromag- 
netic transition is accompanied by a distortion of the 
high-temperature cubic lattice to a monoclinic variant 
with pscudo-rhombohedral geometry. The basic antifer- 
romagnetic structure of MnO can be described in terms 
of a single spin-alignment axis, and its spin dynamics are 
well-understood in terms of simple Heisenberg interac- 
tions. 

We collected neutron total scattering data for a pow- 
dered sample of MnO using the GEM instrument at 
ISIS |l3l |. The experiment spanned a range of momen- 
tum transfers 0.3 < Q < 50 A -1 and was performed at a 
temperature of 100 K. We expected the increased relative 
magnitude of spin displacements at this temperature to 
facilitate our analysis: a study at more "conventional" 
temperatures (e.g., where the spin excitation energies 
are less broad) would place a greater emphasis on res- 
olution effects than proof of concept. Raw data were 
converted to S(Q) and Bragg intensity data for use as 
input for the RMC procedure. We used the RMCProfile 
program to refine an ensemble of ca 600 RMC config- 
urations, each of which contained a 20 x 20 x 20 array 
of fee unit cells (i.e., a 10 x 10 x 10 array of conven- 
tional magnetic unit cells, containing 32 000 Mn atoms). 
Our starting configurations were prepared by assigning 
to each atom a small random displacement from its av- 
erage position; the spin orientations were also varied in a 
similar fashion. For consistency with previous spin-wave 
investigations 0, 0, 0] , we neglected any deviations 
from cubic lattice symmetry; this simplification is known 
to have little effect on the associated analysis. Each pair 
of equilibrium configurations were separated by a mini- 
mum of 150 000 RMC "moves" Q, so that they might 
be considered essentially independent for our analysis. 

A method of extracting spin-dynamical information 
from ensembles of spin configurations will be published 
in detail elsewhere [l5|; there is a strong analogy to the 
known methods of extracting lattice-dynamical informa- 
tion from atomistic configurations [Hi ] . Here we describe 
the most pertinent aspects of the theory. The analysis 
begins by calculating for each configuration (labelled t) 
a number of collective variables of the form 

t(j, k,t) = J||^a + (^,t)exp[ik.r(^)]. (3) 

e 

There will be one of these variables for each of the Z spins 
in the (primitive) magnetic unit cell. Here, Sj is the spin 
quantum number of the spin type j, r(j£) its average 
position in unit cell £, and a + (j£,t) a parameter that 
describes the deviation from the average spin alignment 
axis. All Z collective variables at each wave- vector k are 
assembled into the one column vector i~(k, i). This quan- 
tity is related to the Holstein-Primakoff magnon variable 
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bk(i) [13 y i & the mapping 

V^bk(i) = A(k)-r(M). 



(4) 



Here, the change of basis occurs between the normal 
mode coordinates (the basis of b) and the spin-type 
coordinates (the basis of r), and is given by the (t- 
independent) matrix A, itself constructed from the spin- 
wave mode displacement vectors |l5j . 

We proceed by calculating the i-averaged matrix X, 



E(k) 



: (k)). 



(5) 



This representation is useful because it allows us to ex- 
ploit the orthonormality of the basis for b: in particular, 
the matrix b£(t)b£(f) is diagonal with entries given by 
the mode occupation numbers n(k) ^t|- Substitution of 
Eq. 10} into Eq. (J5J gives 



A(k) • S(k) • A T (k) = h(Kbl) 



(6) 



Since S(b£bJ) is diagonal, its elements [the 7m(k)] are 
given by the eigenvalues ei(k) of E(k). In this way, diag- 
onalisation of S, a matrix that can be constructed (via 
the t) entirely from observed spin displacements in RMC 
configurations, yields the spin- wave occupation numbers. 
The spin-wave frequencies may then be calculated in a 
straightforward manner: 



Wj(k) 



In 



e*(k) 



(7) 



By repeating this analysis for a range of wave-vectors, 
one may construct a set of spin-wave dispersion curves 
from the RMC configurations. The number of k-points 
included in these dispersion curves is limited by the con- 
figurational box size used (which in turn is limited by 
the available computational resources): a box represent- 
ing n a , rib, n c unit cells along axes a, b, c permits the set 
of wave- vectors 

k = ^a* + ^b* + ^c* i a ,i b ,i c eZ. (8) 
n a n b n c 

We proceeded to apply this analysis to our ensemble 
of MnO configurations, calculating spin- wave frequencies 
across the [00£] and [££0] directions of re- 

ciprocal space — these being the axes for which inelastic 
neutron scattering (INS) data were available for compar- 
ison [lfl Hll|. The k-point mesh size (O.lr.l.u.) was lim- 
ited by the number of magnetic — rather than nuclear — 
unit cells in the configuration. Calculation at symmetry- 
equivalent wave- vectors enabled determination of the er- 
rors involved. While we observed two distinct spin- 
wave modes at each wave-vector, inspection of the as- 
sociated mode displacement vectors revealed that their 
relative order in the spin-wave spectrum was not consis- 
tent. This is unsurprising given that the energies of the 
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FIG. 1: (Colour online) MnO spin-wave dispersion curves cal- 
culated from our RMC configurations (data points). Fits 
to the INS data of Peuv lilt (shaded regions, 88.75 K< 
T < 113.75 K) and of Kohgi [nf (red dashed lines, T = 100 K) 
are included for comparison. The relevant directions in re- 
ciprocal space are illustrated on a (110) section of the MnO 
reciprocal lattice. 



two modes are known to be essentially identical at most 
wave- vectors 0, 0, ■ Consequently, the frequencies 
obtained were also averaged over both branches. Our re- 
sults [Fig. ^ are compared here to the dispersion curves 
given in previous INS studies. Pepy gives data for the 
and [OOf] branches at 88.75 K and 113.75 K, 
which serve as upper and lower bounds respectively to 
the expected energy spectrum 0; Kohgi et al report 
data at 100 K, from which they were able to calculate 
dispersion curves along the same three branches together 
with the [££0] direction [ll|. Collins et al. give data at 
only at 4.2 K ^2|> as there is some significant tempera- 
ture dependence in the dispersion, these results are not 
considered further. The curves shown in Fig. ^ are taken 
from the spin-interaction models described in the original 
reports, refined by fitting to the observed INS data. The 
magnon peaks can be relatively broad at such tempera- 
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tures, making absolute energies difficult to assign |lCtlll| . 
Indeed it is interesting that the models differ somewhat 
in their overall energy scale and in the behaviour of the 
spin-wave energies near the zone centre. 

Our purpose in comparing our calculations with the 
results of INS experiments is to establish whether it is at 
all possible to extract a similar level of spin-dynamical 
information from diffraction data. In this sense, it is 
promising that the spin-wave energies we calculate from 
our RMC configurations agree well with the overall en- 
ergy scale observed in INS studies: the level of correlated 
spin "motion" in our RMC configurations is appropri- 
ate for the given temperature and actual spin- wave ener- 
gies. Moreover, there is some clear evidence for variation 
in spin-wave frequency with wave-vector. Certainly, the 
form of the and branches — and to a lesser ex- 
tent that of the [00£] branch — follows closely that of the 
INS-based models. In addition to the general form of the 
dispersion curves, one very encouraging feature is the ap- 
parent increase in spin-wave energy at the M-point rela- 
tive to the zone centre, which reflects the known inequiv- 
alence of parallel and anti-parallel nearest-neighbour cou- 
plings. However, there are regions where the similarity 
between our dispersion curves and those determined us- 
ing INS is less precise: our data along [££0] show rela- 
tively little dispersion, but do occur over an appropriate 
energy range. It is not immediately obvious why the 
diffraction data should be more sensitive to some regions 
of the magnon dispersion than others; what is known 
is that modes at different wave-vectors can contribute 
to differing extents to the real-space distribution func- 
tions 

Despite our employment of state-of-the-art computa- 
tional resources, the RMC configurations prepared were 
not sufficiently large to determine the precise dispersion 
behaviour near the BZ origin. The existence of "gaps" , 
their anisotropy and temperature dependence are often 
important features of spin-wave spectra, and a method 
would have to reflect these to be of general use. However 
a limitation of k-point mesh size is one that can be ex- 
pected to abate as computational power increases. Other, 
more inherent, limitations will exist. For example, it is 
likely that there is some maximum observable spin-wave 
energy: this is certainly true for phonon frequencies 
and magnons suffer additionally from the restricted Q- 
range over which magnetic scattering is observable. 

At this "proof of principle" stage, the important result 
is that it is possible to retrieve even broad features of 
the spin-wave dispersion from powder diffraction data. 
Already, one might argue that relatively sensible param- 
eters for a spin interaction model could be extracted from 
dispersion curves such as those in Fig.^ given that some 
features of the dispersion are well described, and those 
that are less-well fixed by the data might often be con- 
strained by our understanding of spin dynamics. More- 
over, the accuracy of the magnon dispersion curves deter- 



mined from diffraction data can be expected to improve 
significantly. Our methods of treating the very high qual- 
ity data one obtains from instruments such as GEM will 
improve; access to sufficient computational power will in- 
crease. The potential for powder diffraction experiments 
to allow determination of spin- wave excitation spectra is 
itself an important result: it provides a possible mech- 
anism for the exploration of spin dynamics in newly- 
characterised materials, for which single crystal samples 
are not immediately available. 

We acknowledge financial support from EPSRC 
(U.K.), and from Trinity College, Cambridge to A.L.G. 



[1] D. A. Dimitrov, D. Louca, H. Roder, Phys. Rev. B 60, 
6204 (1999); W. Reichardt, L. Pintschovius, Phys. Rev. B 
63, 174302 (2001); M. J. Graf, I.-K. Jeong, D. L. Starr, 
R. H. Heffner, Phys. Rev. B 68, 064305 (2003); I.-K. 
Jeong, R. H. Heffner, M. J. Graf, S. J. L. Billinge, Phys. 
Rev. B 67, 104301 (2003). 

[2] A. L. Goodwin, M. G. Tucker, M. T. Dove, D. A. Keen, 
Phys. Rev. Lett. 93, 075502 (2004); 95, 119901(E) (2005) 

[3] A. L. Goodwin, M. G. Tucker, E. R. Cope, M. T. Dove, 
D. A. Keen, Phys. Rev. B 72, 214304 (2005). 

[4] M. T. Dove, M. G. Tucker, D. A. Keen, Eur. J. Miner. 
14, 331 (2002). 

[5] V. M. Nield, D. A. Keen, Diffuse Neutron Scattering from 
Crystalline Materials (Oxford University Press, Oxford, 
England, 2001); T. Egami, S. J. L. Billinge, Underneath 
the Bragg Peaks: Structural Analysis of Complex Mate- 
rials (Pergamon Press, Oxford, England, 2003). 

[6] R. L. McGreevy, J. Phys.: Condens. Matter 13, R887 
(2001) 

[7] D. A. Keen, M. G. Tucker, M. T. Dove, J. Phys.: Con- 
dens. Matter 17, S15 (2005). 
[8] A. L. Goodwin, M. G. Tucker, M. T. Dove, D. A. Keen, 

Phys. Rev. Lett. 96, 047209 (2006). 
[9] W. L. Roth, Phys. Rev. 110, 1333 (1958); H. Shaked, J. 
Faber, Jr., R. L. Hitterman, Phys. Rev. B 38, 11901 
(1988); D. Hohlwein, J.-U. Hoffmann, R. Schneider, 
Phys. Rev. B 68, 140408(R) (2003). 

[10] G. Pepy, J. Phys. Chem. Solids, 35, 433 (1974). 

[11] M. Kohgi, Y. Ishikawa, I. Harada, K. Motizuki, J. Phys. 
Soc. Jpn. 36, 112 (1974). 

[12] M. F. Collins, V. K. Tondon, W. J. L. Buyers, Int. J. 
Magn. 4, 17 (1973). 

[13] A. C. Hannon, Nucl. Instrum. Meth. A 551, 88 (2005). 

[14] The separation of RMC configurations containing N 
atoms by at least N In 10 "moves" ensures that the po- 
sitions or orientations of at least 90% of atoms or spins 
differ from configuration to configuration. 

[15] A. L. Goodwin, M. T. Do ve, M. G. Tucker, D. A. Keen, 
arXiv:cond-mat/0601559 

[16] M. T. Dove and R. M. Lynden-Bell, Philos. Mag. B 54, 
443 (1986); M. T. Dove, Introduction to Lattice Dynam- 
ics (Cambridge University Press, Cambridge, England, 
1993). 

[17] T. Holstein and H. Primakoff, Phys. Rev. 58, 1098 
(1940). 



